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We discuss formation of an 5-wave bound state in finite volume on the basis of Liischer's phase- 
shift formula. It is found that although a bound-state pole condition is fulfilled only in the infinite 
volume limit, its modification by the finite size corrections is exponentially suppressed by the spatial 
extent L in a finite box L 3 . We also confirm that the appearance of the S- wave bound state 
is accompanied by an abrupt sign change of the S'-wave scattering length even in finite volume 
O ■ through numerical simulations. This distinctive behavior may help us to distinguish the loosely 

bound state from the lowest energy level of the scattering state in finite volume simulations. 
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INTRODUCTION 



00 

In the past few years, several new hadronic resonances have been discovered in various experiments [1]. However, 
, some of the states have unusual properties, which are not well understood from the viewpoint of the conventional 
quark-antiquark or three-quark states. It is a great challenge for lattice QCD to answer the question, whether those 
states are really exotic hadron states. 

We are especially interested in some candidates of hadronic molecular state: the A(1405) resonance as an KN 
£3 I bound state, the /o(980) and ao(980) resonances as S-wave bound states of KK, the X(3872) resonance as a weakly 
"j" 1 ■ bound state of D D , the D s o(2317) and D ,51 (2460) resonances as D^K bound states and so on [2]. Such states lie 
near and below their respective thresholds so that one can view them as "loosely bound states" of two hadrons like a 
deuteron. 

In the infinite volume, the loosely (near-threshold) bound state is well defined since there is no continuum state 
below threshold. However, in a finite box on the lattice, all states have discrete energies. Even worse, the lowest 
energy level of the elastic scattering state appears below threshold in the case if an interaction is attractive between 
r% , two particles [1,0]. Therefore, there is an ambiguity to distinguish between the loosely bound state and the lowest 
■ scattering state in finite volume in this sense. 

Signatures of bound-state formation in finite volume are of main interest in this paper. We may begin with a naive 
question: what is the legitimate definition of the loosely bound state in the quantum mechanics? In the scattering 
theory poles of the S-matrix or the scattering amplitude correspond to bound states. It is also known that the 
appearance of the S'-wave bound state is accompanied by an abrupt sign change of the S'-wave scattering length [f|. 
It is interpreted that formation of one bound-state raises the phase shift at threshold by n. This particular feature 
is generalized as Levinson's theorem Thus, it is interesting to consider how the formation condition of bound 
states is implemented in Liischer's finite size method, which is proposed as a general method for computing low-energy 
scattering phases of two particles in finite volume [3, 0] . 

In this paper, we discuss bound-state formation on the basis of the phase-shift formula in this method and then 
present our proposal for numerical simulations to distinguish the loosely bound state from the lowest scattering state 
in finite volume. To exhibit the validity and efficiency of our proposal, we perform numerical studies of the positronium 
spectroscopy in compact scalar QED model. In the Higgs phase of U(l) gauge dynamics, the photon is massive and 
then massive photons give rise to the short-ranged interparticle force between an electron and a positron exponentially 
damped. In this model, we can control positronium formation in variation with the strength of the interparticle force 
and then explore distinctive signatures of the bound-state formation in finite volume. 
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The organization of our paper is as follows. In SecHH we first give a brief review of Liischer's finite size method 
and discuss bound-state formation on the basis of the phase-shift formula in this method. Sec. IIIII gives details of 
our utilized model, compact scalar QED, and its Monte Carlo simulations. Sees. IIVI and |V]are devoted to discuss 
our numerical results in the 1 5o and 3 Si channels of electron-positron system, respectively. Finally, in Sec. IVI| we 
summarize the present work and give our concluding remark. In addition, there are two appendices. In Appendices 
A, the sensitivity of mass spectra to choice of spatial boundary condition is discussed. We also demonstrate a specific 
volume dependence of the spectral amplitude for either the bound state or the lowest scattering state in Appendix B. 



II. METHODOLOGY 



A. Liischer's finite size method for scattering phase shift 

Let us briefly review Liischer's finite size method [1,0]. So far, several hadron scatteringlengths, eji. tt-tt, tt-K, 
tt-N, K-N . N-N and J/i/i-hadron, have been successfully calculated by using this method jy, 0,11, H, El EH El, El, 

El, EE El E3, El El, M M- 

The total energy of two-particle states in the center-of-mass frame is given by 



E AB (p) = ^m 2 A +p 2 + ^/m 2 B +p 2 , (1) 

where p is the relative momentum of two particles. In a finite box L 3 on the lattice, all momenta are quantized and 
can be labeled by an integer n as p n , which represents the (n + l)-th lowest momentum. Therefore, all two-particle 
states have only discrete energies. 

We introduce the scaled momentum as q = Lp n /2ir with the spatial extent L for periodic boundary condition. 
Although the value of q 2 takes an integer value in the non-interacting case, q 2 is no longer the integer due to the 
presence of the two-particle interaction. This particular feature can be observed through an energy shift relative to 
the energy of the non-interacting two particles, 

AE = E AB (p n )-E AB ( Pn ), (2) 

where the energy of non- interacting two-particle states Eab{Pu) can be evaluated with the quantized momentum p n 
in the free case as p n = 2ir^/n/ L with an integer n. 

It has been shown by Liischer that this energy shift in a finite box with a spatial size L can be translated into the 
S-wave phase shift <5o through the relation [1 [j] : 

TT 3 / 2 J^f 

tan<S (p„) = — ~r at g = Lp n /2n, (3) 
2oo(l,r) 

where the function 2qo(s, q 2 ) is an analytic continuation of the generalized zeta function, Zqq(s, q 2 ) = ^= X)nez 3 ( n2 — 
q 2 )~ s , from the region s > 3/2 to s = 1. The S'-wave scattering length is defined through ao = lim p ^ tan<5o(p)/p. 

If the S'-wave scattering length ao is sufficiently smaller than the spatial size L, one can make a Taylor expansion of 
the phase-shift formula ((3|) around q 2 = 0, and then obtain the asymptotic solution of Eq. ([3|). Under the condition 
2 A and m B , 



p 2 <C m 2 A and m B , the solution is given by 



A w ~ 2lTa ° 



1 + C1 T + C2 Il 



0(L- 6 ), (4) 



which corresponds to the energy shift of the lowest (n = 0) scattering state. The coefficients are C\ = —2.837297 
and C2 = 6.375183 [1, 0|. The reduced mass of two particles /i is given by /i = itia ■ m B /(mA + m B ). An important 
message is received from Eq. ([J} . The lowest energy level of the elastic scattering state appears below threshold on the 
lattice if an interaction is weakly attractive (ao > 0) between two particles. This point makes it difficult to distinguish 
between near-threshold bound states and scattering states on the lattice. 

Here, it is worth noting that the large L expansion formula (j4|) up to 0(L -4 ) gives no real solution of ao for the 
case AE < — ^^hj? EI' w h ue Eq. ([!]) with an expansion up to 0(L~ 4 ) and that up to 0(L~ 5 ) always possesses a 
real and negative solution of ao for AE > 0. A lower bound AE > — 2 \ c ±\iiL' 1 ma y ^ e crucial to identify the observed 
state below threshold as the lowest energy level of the elastic scattering state. 
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For the second lowest (n = 1) scattering state, we also obtain a different asymptotic solution of Eq. ©, which is 
given by a Taylor expansion of the phase-shift formula ([3]) around q 2 — 1 as 

AEq2 = 1 « _ 6ta ^MPl) ^ + c , tan( j o(Sl) + ^ tan 2 Sq{Pi) } + 0(i - 8)j (5) 

where = —0.061367 and c' 2 = —0.354156. Although the sign of tan Jo is not uniquely related to the sign of the 
energy shift, the resulting energy shift AE becomes positive (negative) for the weak repulsive (attractive) interaction 
case (|<5o| ^ 37r/5). Subsequently, one can derive the asymptotic solutions for the higher energy levels of the scattering 
state around q 2 = v > 2 where v = n 2 for integer 3-dim vectors n G Z 3 . For those asymptotic solutions, the 
corresponding relative momentum p n , which we will hereafter abbreviate as p, should vanish as 1/L with increasing 
L. 



B. Bound-state formation in Liischer's formula 



In quantum scattering theory, the formation condition of bound states is implemented as a pole in the S'-matrix or 
scattering amplitude. Therefore, an important question naturally arises as to how bound-state formation is studied 
through Liischer's phase-shift formula ((3]) 



Intuitively, the pole condition of the S'-matrix: S = c 2i& «(p) = ^ot is expressed as 

cot5o(p) = i, 



(6) 



which is satisfied at p 2 = —j 2 where positive real 7 represents the binding momentum. In fact, as we will discuss in 
the following, such a condition is fulfilled only in the infinite volume. However the finite-volume corrections on this 
pole condition are exponentially suppressed by the size of spatial extent L. 

For negative q 2 , an exponentially convergent expression of the zeta function Zqq(s, q 2 ) has been derived in Ref. [22| . 
For s — 1, it is given by 



Z 00 (l,q 2 




-2Try/-q 2 n 2 



(7) 



n£Z 3 



where X)nez 3 means the summation without n = (0,0,0). We now insert Eq. ((7|) into Eq. ^ and then obtain the 



following formula, which is mathematically equivalent to Eq. ([3]) for negative q 2 " 

1 



cot(5 (p) = i + 



1 

2vri 



n£Z 3 



V- 



q 2 n 2 



(8) 



The second term in the r. h. s. of Eq. ijHJ) vanishes in the limit of q 2 



-00. It clearly indicates that negative 



infinite q 2 is responsible for the bound-state formation. Therefore, in this limit, the relative momentum squared p 2 
approaches — 7 2 , which must be non-zero. Meanwhile, the negative infinite q 2 turns out to be the infinite volume 
limit. This representation shows that although the pole condition is fulfilled in the infinite volume, its modification in 
finite volume is described by correction terms, which are exponentially suppressed by the size of spatial extent L oc q. 

Although it was pointed out how the bound-state pole condition could be implemented in his phase-shift formula 
in the original paper another type of large L expansion formula around q 2 = — 00 has been explicitly derived in 
Ref. [2l. 



A£L2__„ = - — 
9 " 00 lii 



12 



1 



7-L 1 - 27(pcot£ )' 



■0(e 



-V2~fL\ 



-0( 7 2 //i 2 ) 



(9) 



where (pcotSo)' — ^p-(pcot<5o)| p 2 = _ 7 2. An ./[/-independent term — 2- corresponds to the binding energy in the infinite 
volume limit. We can learn from Eq. that "loosely bound states" are supposed to receive larger finite volume 
corrections than those of "tightly bound states" since the expansion parameter is scaled by the binding momentum 
7. Furthermore, it can be expected that the bound state of two or more particles has a kinematical nature similar to 
a single particle if the spatial size L is much larger than the size of its compositeness, which may be characterized by 
the inverse of the binding momentum. 



4 



C. Novel view from Levinson's theorem 



At last, a crucial question arises: once the S-wave bound states are formed, what is the fate of the lowest S-wave 
scattering state? The answer to this question might provide a hint to resolve our main issue of how to distinguish 
between "loosely bound states" and scattering states. A naive expectation from Levinson's theorem in quantum 
mechanics is that the energy shift relative to a threshold turns out to be opposite in comparison to the case where 
there is no bound state. Levinson's theorem relates the elastic scattering phase shift Si for the l-th partial wave at 
zero relative momentum to the total number of bound states (Ni) in a beautiful relation [44] ]: 



5,(0)=iV i7 r. 



(10) 



Therefore, if an S-wave bound state is formed in a given channel, the S'-wave scattering phase shift should always 
be positive at low energies. This positiveness of the scattering phase shift is consistent with a consequence of the 
attractive interaction. Conversely, the S'-wave scattering length may become negative (arj < 0) as schematically 
depicted in Fig Q] Consequently, according to Eq. (j4j [45j] , possible negativeness of the scattering length gives rise 
to a positive energy-shift of the lowest scattering state relative to the threshold energy. In other words, the lowest 
(n = 0) scattering state is pulled up into the region above threshold. Therefore, the spectra of the scattering states 
quite resembles the one in the case of the repulsive interaction. If it were true, we can observe a significant difference 
in spectra above the threshold between the two systems: one has at least one bound state (bound system) and the 
other has no bound state (unbound system). 

III. SETUP OF NUMERICAL SIMULATIONS 
A. Compact Scalar QED 

To explore signatures of bound-state formation on the lattice, we consider a bound state (positronium) between an 
electron and a positron in the compact QED with scalar matter: 



(11) 



sites 



which is the compact U(l) gauge theory coupled to both scalar matter (Higgs) fields $ and fermion (electron) fields 
"J. The action of "C/(l) gauge + Higgs" part is described by the compact J7(l)-Higgs model: 



s A M[u,*] = p^[i-si{u g , ltv }]-hY l *i{KVx, lt **+ l *h 

plaq. 



(12) 



link 



where f3 = 1/e 2 and the constraint |$ x | = 1 is imposed. In tree level, the vacuum expectation value of the Higgs field 
and the photon mass are interpreted as (0higgs) ~ a~ x \[h and M pn <~ aT 1 \Jhj ' (3 respectively [24] . In the Higgs phase, 
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one bound state: a <0 




no bound state: a >0 


repulsive: a <0 



FIG. 1: A schematic figure for the scattering phase shift as a function of the relative momentum of two-particle. 
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the Coulomb potential should be screened by the massive photon fields: 

2 p -M ph r 

V(r) ~ . (13) 

4-7T r 

The phase structure of the t/(l)-Higgs model has been well studied on the lattice. Figf3] shows a schematic phase 
diagram of the compact J7(l)-Higgs model. There are three phases: the confinement phase, the Coulomb phase and 
the Higgs phase. The open symbols and filled symbols represent the second-order phase transition points (E: the 
end point {/?, h} = {0.8485(8), 0.5260(9)} {25j and S: the 4-dim XY model phase transition) and the first-order phase 
transition points (T: the triple point {fi,h} ~ {1,0.36} and C: the pure compact U(l) phase transition (3 C ~ 1.01) 
respectively. Lines ET and TC represent the first order line. A dotted line TS corresponds to the Coulomb-Higgs 
transition, of which the order is somewhat controversial in the literature because of large finite size effects. 

B. Monte Carlo simulation 

In this numerical study, we treat the fermion fields in the quenched approximation. Therefore, for update of gauge 
links and Higgs fields, we simply adopt the Metropolis algorithm. First, the acceptance is adjusted to about 30%. 
Then we use 16 hits at each link and Higgs field update. 

Our purpose is to study the S-wave bound state and scattering states through Liischer's finite size method, which 
is only applied to the short-ranged interaction case. Thus, we fix j3 = 2.0 and h — 0.6 in Eq. (TT2"|) to simulate the Higgs 
phase of U(l) gauge dynamics, where massive photons give rise to the short-ranged interparticle force between an 
electron and a positron. We generate U(l) gauge configurations with a parameter set, (/?, h) — (2.0,0.6), on L 3 x 32 
lattices with several spatial sizes, L = 12, 16, 20, 24, 28 and 32. Statistics for each volume calculation are summarized 
in Tabled] 



TABLE I: Simulation statistics 



Spatial size (L) 


12 


16 


20 


24 


28 


32 


# of conf. 


960 


1920 


1280 


720 


720 


480 



Once the parameters of the compact Z7(l)-Higgs action, (/3,h), are fixed, the strength of an interparticle force 
between electrons should be frozen on given gauge configurations. However, if we consider the fictitious Q-charged 
electron, the interparticle force can be controlled by this charge Q since the interparticle force is proportional to 
(charge Q) 2 . Within the quenched approximation, this trick of the Q-charged electron is easily implemented by 
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FIG. 2: Schematic phase diagram of the compact [/(l)-Higgs model in the fixed modulus case. A star mark represents our 
simulation point as (/3,h) — (2.0, 0.6). 
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replacing U(l) link fields as 



into the Wilson-Dirac matrix: 



X.fl 



D w [U Q Uy = S x , y - k£ [(1 - 7/*)0*+/*.» + (1 + 



x-fi,y 



(14) 



(15) 



15 



where k is the hopping parameter. 

For the matrix inversion, we use the BiCGStab algorithm [26[ and adopt the convergence condition \r\ < 10 
for the residues. We calculate the electron propagators (0\^ (x)^ (y)\0) ~ D^\U^] XiV with both periodic and anti- 
periodic boundary conditions in the temporal direction. Then, we adopt the avera ged p ropagator over the boundary 
conditions. This procedure provides an electron propagator with 2T-periodicity [271 l28j|. 



C. Spectrum of single electron 



To evaluate a threshold energy of the electron-positron (e e + ) system, it is necessary to calculate the electron mass 
nonperturbatively by the following two-point correlator, 



G e (*;Pn) = Z¥ ^Tr{P + <0|vl/(x,i)vI/(y,0)|0)e i P-( x -^}, (16) 
x,y 

where V+ = 1H 2 74 and p„ = with neZ 3 for the periodic boundary condition in spatial directions. Here, we have 
set the lattice spacing to unity (a = 1). This electron two-point correlator is gauge-variant, so gauge fixing is required. 
We fix to the Landau gauge. However, it is well known that the pure compact U{\) gauge theory in the Coulomb 
phase leads to a serious problem of the Gribov ambiguity in the gauge-fixing procedure. We adopt the modified 
iterative Landau gauge fixing, which is proposed in Ref . [29j , to avoid the Gribov copy effect on gauge-variant electron 
correlators as much as possible. Here, we remark that the Gribov ambiguity is not observed to be severe in the Higgs 
phase of compact scalar QED, where our simulations are performed, as is also true in the confined phase [29| . 



TABLE II: Two parameter sets (Q, k) for electron fields and resulting rest masses of a single electron in lattice units. 



charge Q 


K 




3 


0.1639 


0.479036(75) 


4 


0.2222 


0.50396(59) 



1. Volume dependence of electron rest mass 

According to our previous pilot study (30l | , numerical simulations are performed with two parameter sets for fermion 
(electron) fields, (Q, k)=(3, 0.1639) and (4, 0.2222), which are adjusted to yield almost the same electron masses 
M e ss 0.5 for both charges. First, we calculate the electron mass at rest (po = (0, 0, 0)). The electron mass is obtained 
by a single exponential fit, which takes into account the 2T-periodicity in our simulations, to the two-point correlator 
of a single electron (fTB"|) . In Fig[3j we show the volume dependence of the electron mass for three-charged (the left 
panel) and four-charged (the right panel) electron fields. In both cases of Q = 3 and 4, there is no appreciable finite 
size effect on the electron mass if the spatial lattice size L is larger than 16. We take a weighted average of the five 
masses in the range 16 < L < 32 to evaluate values in the infinite volume limit, which are hereafter used in estimating 
a threshold energy of two-electron states. In FigJ31 solid horizontal lines represent the average values taken as the 
infinite volume limit, together with their one standard deviation (dashed lines). A summary of the infinite- volume 
values is given in Table HU 
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Q=4 



10 15 20 25 30 35 10 15 20 25 30 35 

L L 



FIG. 3: Measured electron masses as a function of spatial size L in lattice units. The left (right) panel is for three-charged 
(four-charged) electron at k = 0.1639 (0.2222). The horizontal solid line in each panel represent the value of M e in the infinite 
volume limit, which is obtained from the weighted average of five data, with its la deviation (dashed lines). 



TABLE III: Fitted masses of single electrons (Q — 3,4) with zero momentum and non-zero lowest momentum (pi = 2n/L) at 
six different lattice volumes L 3 x 32. 



Spatial size 


Q = 3 




Q = 4 




L 


M e 


E.(px) 




Ee(pi) 


12 


0.48091(53) 


0.67364(85) 


0.5050(38) 


0.6941(55) 


16 


0.47916(21) 


0.60248(34) 


0.5031(14) 


0.6235(18) 


20 


0.47889(18) 


0.56255(26) 


0.5057(12) 


0.5861(13) 


24 


0.47916(18) 


0.53953(26) 


0.5029(12) 


0.5610(13) 


28 


0.47892(15) 


0.52485(24) 


0.5041(14) 


0.5483(14) 


32 


0.47912(15) 


0.51506(27) 


0.5036(14) 


0.5388(16) 


oo 


0.479036(75) 




0.50396(59) 





2. Dispersion relation 



Next, we examine the dispersion relation of the single electron in our simulations in order to study the effects of the 
finite lattice spacing. We calculate the electron correlation (fl"6|) with non-zero lowest momentum, pi = ^(1,0,0), to 
measure the energy level of the non-zero momentum single electron. All measured values are tabulated in Table ITLT1 
In FigHl we compare our measured energies E e (pi) at several spatial lattices with a couple of theoretical curves, which 
are evaluated from two types of the dispersion relation with the measured rest mass: the continuum-type dispersion 
relation 



K on (Pn) = V^ + pI (17) 

and the lattice dispersion relation for free Wilson fermions [3l[ 



E^iPn) = COsh- 1 + (l->/^g + y) 2 +Pn f (18) 

where W — e Me — 1, p„ = ^(n x , n y ,n z ), and pj; = X)fe sm2 [T" nfc ]- The son d curves obtained from the lattice 
dispersion relation are clearly closer to the measured energies in both Q = 3 and Q = 4 cases. The finite lattice 
spacing effects on the single electron spectra are not negligible even at the lowest momentum. Recall that the relative 
momentum of two particles is a key ingredient when we determine the scattering phase shift from Eq.Q. In this sense, 
the lattice dispersion relation is preferable so as to reduce the systematic error stemming from the lattice spacing 
artifact in determination of the relative momentum of two-particle states. Through out this paper, we use the lattice 
dispersion relation (|18p in the analysis of the scattering phase shift through Liischer's formula ([3]). 



s 




FIG. 4: The electron energy at the non-zero lowest momentum (pi = 2n/L) as a function of the spatial size L in lattice units. 
The left (right) panel for the Q — 3 (Q = 4) case. Full circles represent the measured values. The solid (dashed) curves are 
theoretical curves evaluated from the lattice (continuum-type) dispersion relation with the measured rest mass. 



D. Diagonalization method 

We are especially interested in the 1 S'o and 3 Si states of the e~e + system, where the electron-positron bound 
state (positronium) could be formed even in the Higgs phase, and 3 Si positronium are described by the bilinear 
pseudo-scalar operator ^ x 75^x and vector operator ^> x l^x respectively. Therefore, we may construct the four-point 
functions of electron-positron states based on the above operators. We arc interested in not only the lowest level of 
two-particle spectra, but also the 2nd and 3rd lowest levels. In order to extract a few low-lying energy levels of 
two-particle system, we utilize the diagonalization method proposed by Liischer and Wolff [32|]. We consider three 
types of operators for this purpose: 

X 

n w(t) = ^£*(y>*) r *( x '*)' (20) 

^M(t) = -^^(y^r^x,^ 1 '^, (21) 

where pi = ^(1,0,0) and T — 75 (7^) for the 1 1 5 , ( 3 5i) e _ e + state. The first operator corresponds to a simple 
local-type operator where only the total momentum of two particles is fixed to be zero, but both the electron and 
the positron can carry non-zero relative momentum under the total momentum conservation. The second operator 
projects both the electron and the positron onto zero momentum, while the relative momentum of the e~e + system 
is constrained to the non-zero lowest momentum (pi = |pi| = ^7) in the third operator. Therefore, we can expect 
that each type of operators has better overlap to a specific two-particle state: n = and n = 1 scattering states have 
strong overlap with fl\y and respectively, while the bound state has the better overlap with flp than Qw and 

We construct the 3x3 matrix correlator from above three operators 

G y (t) = <0|^(*)fit(0)|0) (22) 
and then employ a diagonalization of a transfer matrix M(t,to), which is defined by 

M(t, t ) = G(t )- 1/2 G(t)G(t )" 1/2 , (23) 

where to is a reference time-slice. If only three states are propagating in the region t > to, the energies of three 
two-particle states E a (E2 > -Ei > Eq) are given by the eigenvalues of M(t, to): 

X a (t,to)^e-^- t0 ^ (a = 0,1, 2), (24) 

where E a is independent of to- An assumption that three low- lying states become effectively dominant for an ap- 
propriately large time-slice to, can be determined by checking the sensitivity of E a with respect to variation of the 
reference time-slice to- 
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In this study, the random noise method is employed to calculate i7 p source operators in Eq. with the number of 
noises taken to be one. Technical details of this method are described in Ref. [ill. H3|. We note that all contributions 
from disconnected diagrams in Eq. (|22j) are simply ignored in our numerical calculations. 

IV. NUMERICAL RESULTS IN THE CHANNEL 

In this section, we focus on numerical results in the 1 Sq channel of the e~e + system. Results obtained in the 3 Si 
channel will be separately discussed in the next section. 

A. Ground state of So 

Let us begin with the ground state in the channel. It is not necessary to employ the diagonalization method 
for the spectroscopy of the ground state. We first show the effective mass plot for two diagonal components of the 
3x3 matrix correlator. Figs. [5] show the effective mass of the PP correlator and the WW correlator in simulations 
at spatial extent L = 28 for Q = 3 (left panel) and Q = 4 (right panel). At a glance, there are apparent operator 
dependencies. A very clear plateau appears for the WW correlator in the Q = 3 case, while the same quality shows 
up for the PP correlator in the Q — 4 case. This drastic change in operator dependence strongly suggests a signature 
of bound-state formation in the Q = 4 case, since the WW correlator is expected to have a large overlap with the 
lowest (n = 0) scattering state rather than the bound state. In addition, the energy of the 1 5'o ground state is close 
to the threshold energy (2M e ~ 0.958) in the case of Q = 3, while there is a large energy gap between the ground 
state energy and the threshold energy (2M e ~ 1.008) in the case of Q = 4. Therefore, we may naively conclude that 
the ground state in Q = 4 is the positronium state . 

To make a firm conclusion on this point, we next show the volume dependence of the ground state energy in Fig. [5] 
In the left panel (Q — 3), we plot ground state energies measured at each L together with the threshold energy as 
horizontal lines, which are estimated by 2M e and its one standard deviation. An upward tendency of the L dependence 
toward the threshold energy is clearly observed as spatial size L increases. We also include a lower bound for the 
asymptotic solution of the scattering state. All data points are located well above this lower bound. From those 
observations, we can conclude that the observed ground state in the Q = 3 is definitely the lowest (n = 0) scattering 
state. 

On the other hand, in the right panel (Q = 4), all data points are located far below the threshold energy and also 
the lower bound for the asymptotic solution of the scattering state. Indeed, data are well fitted by the form: 

E(L)=A+jeM-lL), (25) 

which is inspired by the asymptotic solution of the bound state, Eq. ©. Finding A ^ 2M e directly indicates 
that the energy gap from the threshold remains finite in the infinite volume limit. We perform two types of fitting 
procedure with this form. First, a full three-parameter fit is employed. Second, we take into account a relation 
between two parameters A and 7 according to Eq. The parameter A is the value of ground-state energy in 

the infinite volume, while 7 corresponds to the binding momentum related to the pole location of the ^-matrix as 
p 2 = — 7 2 (< 0). Therefore, an explicit constraint between two parameters A and 7 can be imposed through the 
relation 7 = \J M 2 — A 2 /A referred to the measured electron mass M e . Then a two-parameters fit is carried out. All 
fitting results are tabulated in Table. ITVl Either procedure provides reasonable fits with about x 2 /d.o.f. ~ 1. The 
resulting values of A in both fits are approximately consistent with each other, while some differences show up in 
other parameters. 

We here stress that the obtained value of A is significantly far from the threshold value 2M e ~ 1.008 and therefore 
the energy gap AE = A — 2M e clearly remains finite in the infinite volume limit. A bound state of an electron and 
positron is certainly formed in simulations with charge-four electrons. 



TABLE IV: Fitting results for the So ground state of Q = 4 electron fields using a fitting form of Eq. (|25[) . In the table, "Fit 
1" and "Fit 2" stand for the fully three-parameters fit and the two-parameters fit with a constraint between A and 7. 





A 


B 


7 


X 2 /d.o.f. 


Fit 1 


0.76395(26) 


-0.78(30) 


0.068(28) 


0.87 


Fit 2 


0.76350(11) 


-3.92(89) 


constrained as WMi — A 2 / 4 


1.58 
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FIG. 5: The effective masses in the So channel on the lattice with L = 28 as a function of the time-slice t in lattice units. 
The left (right) panel is for Q — 3 (Q = 4) electron fields. Full circles (full squares) symbols are obtained from the PP (WW) 
correlator. 
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FIG. 6: Energies of ground states in the 1 So channel of the e~e + system as a function of spatial lattice size L. In the left panel 
(Q — 3), full circles are measured energies at each lattice size. Horizontal lines represent the threshold energy 2M e and its 1 
standard deviation (dashed lines). The dashed curve shows the lower boundary for the convergence of the large L expansion 
formula around q 2 = as Eq. @. In the right (Q = 4) panel, full circles are measured energies at each lattice size, and the 
solid (dashed) curve is the fitting result by using a form E = A + B exp(— yL)/L, which is inspired by the large L expansion 
formula at infinite negative q 2 as Eq. @, with (without) a constraint between two parameters (A and 7). 



B. Excited state of 1 Sq 



In the previous subsection, we confirm that the simulations with three-charge electrons provide the purely elastic 
scattering system without bound states (unbound system), while four-charge electrons give rise to at least one bound 
state as the ground state in the 1 5'o channel of the e~e + system (bound system). We now can explore the difference 
of spectra between the unbound system (Q = 3) and the bound system (Q = 4). 

We calculate the eigenvalues X a (t,to) of the transfer matrix M(t,to) for to = 7 at all L except for L = 32 where 
to = 9 is chosen. First we show the effective mass plots for all three eigenvalues X a (t, to = 7) in simulations at L = 28 
in Fig. [71 The diagonalization method with our chosen three operators successfully separates the first excited state 
and the second excited state from the ground state. 

In the left panel (Q = 3), the lowest and the second lowest states show very clear plateaus started from t ~ 5, 
which is earlier than our reference time-slice tg- The ground state and the first excited state correspond to the lowest 
(n = 0) scattering state and the second lowest (n — 1) scattering state. Those two-particle energies E® e and El e 
are close to twice the single electron energies, 2E e (p ) = 2M e and 2E e (p 1 ), respectively. Needless to say, the energy 
of the lowest state in the diagonalization method is consistent with the energy obtained by the WW correlator. By 
detail analysis of the spectral amplitude (see, Table fVl and Appendix B), we confirm that the WW correlator and 
the MM correlator are dominant in Ao and Ai respectively as expected. Although the effective mass of the third 
eigenvalue A2 gradually approaches some plateau around t « 20, statistical errors becomes large in the plateau region. 
A2 is dominated by the PP correlator, which can overlap with any relative momentum scattering state, so that the 
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contamination from higher relative momentum (n > 3) scattering states is inevitable in the earlier time-slice. 

For the bound system (Q — 4) , all three eigenvalues show clear plateaus started from t w to = 7 in the effective mass 
plot. Again, the energy of the lowest state in the diagonalization method agrees well with the one obtained from the 
PP correlator. The obtained eigenvectors also indicate that the PP correlator is dominant in the Ao eigenvalue, while 
the second and third eigenvalues are mostly composed of the WW correlator and the MM correlator, respectively. 
As we mentioned, the PP correlator possibly has overlap with any relative momentum scattering states. However, 
here, the PP correlator has dominant overlap with the bound state as shown in Table |V] This is because the spectral 
weight of two-particle states relative to the single particle state, such as a bound state, could be suppressed in the 
PP correlator by an inverse factor of the volume, 1 /L 3 [46[ . 

Finally, a summary table of low- lying spectra in the 1 Sq channel in simulations at L = 28 is given in Table IVll 



TABLE V: Summary of the normalized spectral weights (A a )i in the So channel on the lattice with L = 28. A definition of 
the normalized spectral weights is described in Appendix B. 







ground state 


1st excited state 


2nd excited state 


charge Q 


operator i 


a = 


a = 1 


a = 2 


3 


P 


0.02166(72) 


0.1324(96) 


0.846(10) 




W 


0.99735(37) 


0.00177(11) 


0.00088(31) 




M 


0.001155(48) 


0.9907(11) 


0.0081(11) 


4 


P 


0.9597(83) 


0.0064(21) 


0.0339(81) 




W 


0.01567(69) 


0.9841(14) 


0.00018(78) 




M 


0.0521(15) 


0.0023(17) 


0.9456(24) 





TABLE VI: Summary of low-lying spectra in the 


So channel on the lattice with L 


= 28. 




charge Q 


eigenstate a 


energy E a 




kinds of state 


3 


ground state (a = 0) 
1st excited state (a = 1) 


£° e =0.95672(31) 
£^=1.04062(41) 


n — 
n — 


scattering state 

1 scattering state 


4 


ground state (a = 0) 
1st excited state (a = 1) 
2nd excited state (a = 2) 


M bs = 0.76370(17) 
E° e = 1.0119(26) 
Elz= 1.1044(53) 


n — 
n — 


bound state 

scattering state 

1 scattering state 



• ground 

■ 1st excited 

♦ 2nd excited 



Q=3 L=28 



» ground 
■ 1 st excited 
► 2nd excited 



A; 



Q=4 L=28 



FIG. 7: The effective mass plots for each eigenvalue X a (t, to) of the transfer matrix defined in Eq. (J23J at a reference time-slice 
to — 7 on the lattice with L = 28. Full circles, squares and diamonds represent the ground state (a = 0), the first excited state 
(a = 1) and the second excited state (a = 2). The left (right) panel is for Q = 3 (Q — 4). 
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C. Distinctive signatures of bound-state formation 

1 . Sign of energy shift 

Suppose that Luscher's finite size method reflects all of the essential nature of the scattering theory in the quantum 
mechanics; formation of the S'-wave bound state is accompanied by an abrupt sign change of the scattering length. 
Thus, we can expect that the second lowest energy state, which corresponds to the lowest (n = 0) scattering state, 
should be located near and above the threshold energy (2M e ) if a bound state is formed. This is quite in contrast 
with the case if there is no bound state: the second lowest energy state, which should be the n = 1 scattering state, 
is located near below (above) the energy level of non-interacting two-particle system with non-zero lowest momentum 
as 2E e (pi) in the attractive (repulsive) channel. 

Here we show our observed L-dependence of the energy level of the second lowest state in Fig. [8] The data plotted 
appear in Table IVII1 In the left panel (Q = 3) , measured energy levels are very close to the n = 1 threshold energy 
which is given by twice the single electron energy at non-zero lowest momentum p\. As the spatial size L increases, 
the energy levels approach this n = 1 threshold energy from below. This is consistent with a behavior of the n = 1 
scattering state predicted by Eq. (J3J) for the weakly attractive interaction without bound states. Therefore, one can 
identify the second lowest energy state as the n = 1 scattering state for Q = 3. 

In the right panel (Q = 4), an expected feature comes out. The horizontal line represents the n = threshold 
energy estimated by twice the electron rest mass. Clearly, the energy levels of the second lowest state approach this 
n = threshold energy from above. The energy shift from the threshold vanishes as the spatial size L increases. 
Therefore, the second lowest energy state must be the n = scattering state. It is worth emphasizing that the sign of 
AE = E® e — 2M e is opposite in the case of Q = 3 where there is no bound state. Of course, this sign is directly related 
to the sign of the S-wave scattering length. Thus, our numerical simulations show that formation of the S-wave bound 
state is really accompanied by an abrupt sign change of the scattering length. 

Furthermore, in Fig. [51 the volume dependence of the energy level of the third lowest state in the Q = 4 case 
shows the "repulsive" feature as the n = 1 scattering state even in the attractive channel. This is attributed to the 
consequence of Levinson's theorem, which allows the case, tan Jo < 0, for the attractive interaction. 

What is surprising here is that one of the most important features, namely Levinson's theorem, in the quantum 
scattering theory is inherited in Luscher's finite size formula. Meanwhile, we realize what is a proper signature of 
bound state formation in finite volume on the lattice. Even in a single simulation at fixed L, we can distinguish the 
near-threshold bound state from the lowest (n = 0) scattering state through determination of whether the second 
lowest state appears just above the threshold or near the n = 1 energy level of non-interacting two-particle states. 



TABLE VII: Energies of the ground state, the first excited state and the second excited state (only for Q = 4), which are 
obtained by a single cosh fit with 2T periodicity, in the 1 So channel of the e~e + system at five different lattice volumes L 3 x 32. 



Spatial size 


Q = 3 




Q = 4 






L 




E\e 


M hs 


Eee 




16 


0.95183(40) 


1.16903(55) 


0.76231(26) 


1.0242(52) 


1.287(33) 


20 


0.95445(36) 


1.10435(36) 


0.76300(22) 


1.0145(40) 


1.197(14) 


24 


0.95644(37) 


1.06590(36) 


0.76314(22) 


1.0126(34) 


1.1318(93) 


28 


0.95672(31) 


1.04062(41) 


0.76370(17) 


1.0119(26) 


1.1044(53) 


32 


0.95767(32) 


1.02440(49) 


0.76359(20) 


1.0113(32) 


1.0888(68) 



2. Bound-state pole condition 

As we discussed in Sec. Ill B[ the formation condition of the S'-wave bound state, cot Jo(p) = i, is definitely im- 
plemented in Luscher's phase-shift formula ([3]) at negative infinite q 2 , which corresponds to the limit of L — > oo. 
According to the original paper pj], for negative q 2 , we introduce the phase cto(k), which is defined by an analytic 
continuation of Sq into the complex p plane through the relation 

tant7o( K ) = — i tan Jo (p), (26) 

where k — —ip. Therefore, the bound-state pole condition in the infinite volume reads cot 00(7) = ~ 1 f° r the binding 
momentum 7 Q. As we pointed out in Sec. IIIB1 the finite volume correction on this condition is exponentially 
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FIG. 8: Energies of the first excited state in the 1 So channel of the e~e + system as a function of spatial lattice size L. The 
left (right) panel is for Q — 3 (Q — 4). The horizontal line represents the threshold energy determined by 2M e . Full circles 
are measured energies for the first excited state. Solid curves with full squares shows twice of the single electron energy with 
non-zero smallest momentum pi = 2n/L. 




FIG. 9: Energies of the second excited state in the 1 S'o channel of the e e" 1 
Q = 4. All symbols are defined as in Fig. [8] 



system as a function of spatial lattice size L for 



suppressed by the spatial extent L in a finite box L 3 : 



lim cot ao(ft) 

K^7 



N, 



-1 + — le 



-L 7 



0(e 



-V2L-, 



(27) 
(28) 



where the factor N„ is the number of integer vectors n £ Z with v = n 2 . Therefore, if the bound state is formed, we 
may observe the phase cxq satisfies lim K ^ 7 ctq(k) — — j — e L where e^(> 0) vanishes as the spatial size L increases. 

We want to examine this bound-state pole condition numerically in the known bound system. As described pre- 
viously, it is found that our simulation in the Q = 4 case yields an S'-wave bound state as the ground state in the 
1 S'o channel. Thus, we determine the phase <j from an energy level of the ground state in the Q = 4 simulation 
by using Liischer's formula ((3]). We first calculate the relative momentum of two particles (electron-positron) from 
the measured energy level of the ground state E ee by matching with twice the single electron energy 2E e (p). As we 
discussed in Sec lIII C"2l we prefer to use the lattice dispersion relation (fT8|) for a formula of the single electron energy 
E e (p) in order to avoid lattice discretization errors as much as possible. 

In Fig. [TD1 we plot the phase cto in the 1 5'o channel as a function of p 2 (< 0). The data plotted appear in Table IVLLTl 
One can easily observe that the phase <tq approaches —45 deg. (— 7r/4) from below as the spatial size L increases. 
Even at the smallest spatial extent L — 16, the phase ao is very close to — ir/4. Needless to say, observed values of 
p 2 (= — 7 2 ), which are related to the binding energy of the bound state, are almost insensitive to the spatial size L 
within statistical errors. Thus, we confirm that our observed "bound state" in finite volume approximately fulfills the 
pole condition of the S'-matrix. 
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A more rigorous way to test for bound-state formation would be to use an asymptotic formula for the finite volume 
correction to the pole condition as Eq. (|28|) . In Fig. [Til we plot the value of cotco versus the spatial lattice extent L 
and show two fit results using Eq. ((28]) with different numbers of exponential terms (one term and three terms). As 
shown in Table \TK\ the optimum number of exponential terms, which yields a convergent result of 7, is about three. 
However, results with one term and three terms are quite consistent with each other because of rapid convergence. 
Then both fit curves in Fig. [TTJ reproduce all data points very well. 



TABLE VIII: Summary of the relative momentum squared p 2 , the phase ao and cot ao measured in the 1 So channel for Q — 4 
at five different lattice volumes L 3 x 32. 





16 


20 


24 


28 


32 


P 2 

a (deg.) 

COt (Tq 


-0.1252(17) 
-45.1218(58) 
-0.99576(20) 


-0.1281(15) 
-45.02012(99) 
-0.999298(35) 


-0.1245(15) 
-45.00444(26) 
-0.9998449(91) 


-0.1259(17) 
-45.000860(65) 
-0.9999700(23) 


-0.1252(18) 
-45.000186(17) 
-0.99999352(58) 



D. e e + elastic scattering phase shifts 

Finally, we evaluate the elastic scattering phase shift of both the unbound system (Q = 3) and the bound system 
(Q = 4) using Luscher's formula ([3]). 
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▲ L=28 




L=32 
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FIG. 10: Phase ao as a function of relative momentum squared p 2 in the Q = 4 simulation. Different symbols represent the 
values obtained from simulations with different spatial lattice sizes. 
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FIG. 11: cot ao in the So channel for Q = 4 as a function of the spatial lattice size L. The solid (dashed) curve represents a 
fitting result using Eq. (|28[l with only a leading exponential term (three exponential terms). 
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TABLE IX: Fitting results for the bound-state pole condition in the 1 So channel for Q — 4 using Eq. (|28[) with variation of 
the number of exponential terms. 



fitting range (L) 


# of exp. terms 


7 


X 2 /d.o.f. 


16-32 


1 


0.3524(11) 


1.79 




2 


0.3547(11) 


0.86 




3 


0.3549(11) 


0.85 




4 


0.3549(11) 


0.85 



1. Unbound system (Q = 3) 

In the Q = 3 case, as we described previously, there is no bound state. The ground state and the first excited 
state correspond to n = and n = 1 scattering states respectively and those energy levels are successfully separated 
by the diagonalization method. Then we can measure the scattering phase shifts So(p) at two different kinematical 
points, which correspond to the relative momenta of the two particles (electron-positron) p for both n — and n = 1 
scattering states. However, as for the n = scattering state, the relative momentum squared is negative (p 2 < 0) 
because of the attractive interaction between the electron and the positron. Therefore, we only access the phase Co 
from the energy level of the n = scattering state in this sense. However, we consider the effective-range expansion 
for the scattering phase as pcot 5o(p) = + \ r oP 2 + 0(p 4 ) in the vicinity of zero relative momentum. We then 
assume this expansion is still valid for negative p 2 . Therefore, we may translate the phase <jq to the scattering phase 
shift in the following relation 

lim KCotero(K) = limp cot 5o(p) = — . (29) 

In other words, we approximately identify the value of ctq(k) at p 2 = —k 2 to the scattering phase shift 5q(p) at 
p 2 = +k 2 . On the other hand, the relative momentum squared of the n = 1 scattering state is definitely positive so 
that we directly access the scattering phase shift So through Liischer's formula without any approximation [47| . 

We plot scattering phase shifts measured from both n = and n = 1 scattering-state energies in the Q = 3 case as 
a function of the relative momentum squared in the left panel of Figs. Q2] and in the left panel of Figs. [13] following 
the above descriptions. 



2. Bound system (Q — 4) 

In the Q = 4 case, it is found that the ground state corresponds to the 5-wave bound state. As we discussed in 
Sec. IIII B| the first and the second excited states should be n — and n = 1 scattering states, respectively. Although 
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FIG. 12: Scattering phase shifts So from the energy level of the n — scattering states (corresponding to the ground state for 
Q — 3 and the first excited state for Q — 4). The horizontal axis is the squared relative momentum of the e~e + system. The 
left (right) panel is for Q — 3 (Q — 4). Different symbols represent the values obtained from simulations with different spatial 
lattice sizes. Remark that the sign of So for Q = 4 is opposite to that for Q = 3. 
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FIG. 13: Scattering phase shifts So from the energy level of the n = 1 scattering states (corresponding to the first excited state 
for Q — 3 and the second excited state for Q = 4). The horizontal axis is the squared relative momentum of the e~e + system. 
The left (right) panel is for Q = 3 (Q = 4). Different symbols represent the values obtained from simulations with different 
spatial lattice sizes. Remark that the sign of So for Q — 4 is opposite to that for Q — 3. 



we can not access any information of the scattering phase shift from the energy level of the ground state, we instead 
determine the scattering phase shift from the energy levels of the first and second excited states. In contrast to the 
purely scattering system without bound states {Q = 3), the relative momentum squared is given as the positive value 
even from the lowest (n = 0) scattering state, which appears above the threshold. Therefore, we can simply evaluate 
the phase shift Sq using Liischer's formula with measured p 2 . 

In the right panels of both Figs. [T2] and Figs. [T31 we plot our measured scattering phase shifts for Q = 4 ver- 
sus the relative momentum squared. Here, the values of the phase shift Sq are simply restricted to the interval 
(— ^ , -|]. Therefore, we observe negative phase shift Sq despite an attractive interaction between the electron and the 
positron. Roughly speaking, the phase shift So monotonically increases as p 2 decreases and approaches zero toward 
zero momentum squared. It implies that the S-wetve scattering length ao is negative for the bound system (Q = 4). 



3. Scattering amplitude 

Here, we define an S*-wave scattering amplitude T(p) 

T( P ),^M^, (30) 
P 2 

where E ee represents the measured energy of the scattering state. Analyticity of the scattering amplitude T(p) allows 
us to consider the following fit ansatz: 

T{p) = d + d lP 2 + d 2 p i + d 3P 6 + d iP s , (31) 

which is a simple polynomial function in the relative momentum squared p 2 . The results of the fit are summarized in 
Table. IXl For Q = 4, a linear fit with respect to p 2 is enough to describe the data with reasonable x 2 /d.o.f., while the 
fourth order polynomial fit still yields large x 2 /d.o.f. in the case of Q = 3. The latter point will be discussed before 
this session is closed. 

We then obtain a global p 2 -dependence of the phase shift Sq in the measured region of p 2 , which is deduced from 
fitted results of the scattering amplitude. We show all measured phase shifts Sq, which are obtained from the energy 
levels of both n = and n = 1 scattering states, in Figs. Q~J] Solid curves represent inferred p 2 -dependence of the 
phase shift with the band of their errors. 

In the right panel (Q = 4), we take into account the modulo- n ambiguity in determination of the phase shift <5o 
because of the bound system and raise the phase shift by an additional 7r in order to fulfill Levinson's theorem. That 
is why the phase shift data starts from it and monotonically decreases as p 2 increases. All data are well covered with 
rather wide bands of error associated with the global fit. 

On the other hand, in the left panel (Q = 3), two data sets determined from energy levels of n — and n = 1 
scattering states seem not to be smoothly connected with each other due to the lower data points from the n = 1 
scattering state at L = 28 and 32. We remark that although statistical errors on all points are rather small, a hidden 
and large systematic error stems from an order Ola) lattice artifact in the determination of p 2 . As we discussed in 
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Sec. IIII C"2l we have used the lattice dispersion relation in the analysis of the scattering phase shift. The continuum- 
type dispersion relation yields smaller estimations of p 2 than those obtained from the lattice dispersion relation. These 
differences are far beyond statistical errors, especially for p 2 obtained from the n = 1 energy level in the Q = 3 case. 
Furthermore, discrepancies are largely enhanced in determination of the scattering phase shift through the Luscher 
finite formula. The scattering phase shift from the n = 1 energy level for Q — 3 typically increases by about a factor 
of two, if the continuum-type dispersion relation is utilized in the whole analysis. 
At the low-energy limit, the scattering amplitude becomes 

lim T(p) = a Q M e . (32) 

Therefore, the fitting parameter do m Eq. (|3"Tj) is associated with the scattering length clq. We then obtain the 
scattering lengths as ao = 1.46(5) for Q = 3 and —2.28(40) for Q = 4 in lattice units, which are much smaller than 
our utilized lattice sizes (L > 12). Needless to say, the sign of the scattering length for Q — 4 is opposite to that for 
Q = 3 due to formation of one bound state in the case of Q — 4. 



TABLE X: Fitting results for the scattering amplitude for Q = 3 and Q — 4 using Eq. (|31[) . 



charge Q 


do 


di 


d 2 


d 3 


di 


X 2 /d.o.f. 


3 


0.697(26) 


-33.8(7.4) 


102(25) x 10 


-103(28) xlO 2 


34(10) xlO 3 


7.26 


4 


-1.15(20) 


-4.6(7.3) 








0.52 




FIG. 14: All measured scattering phase shift So as a function of the squared relative momentum. As for the Q = 4 case (the 
left panel) where the formation of bound state is observed, the scattering phase shifts 80 are shifted as So — ► 80 + 7r according 
to Levinson's theorem. Solid curves represent the fitting results with the band of their errors. 



V. NUMERICAL RESULTS IN THE 3 Si CHANNEL 
A. Low-lying spectra in the Q — 4 case 

The 3 Si spectroscopy has been done in exactly the same way as the So case by using the bilinear vector operator 
^xl^x- As for the Lorentz indices, we take an average over the spatial indices so as to gain possible reduction of 
statistical errors. After we perform the diagonalization method with the 3x3 matrix correlator constructed with 
three operators in Eq. (j2"Tj) . we get the energy spectra of both the ground state and the first excited state. 

In the Q = 4 case, we have concluded that one bound state is formed in the 1 5'o channel as described in the previous 
section. The binding energy B = \M}> S — 2M e \ is rather large as B ss M e /2. The observed bound state should be 
a "tightly bound state" rather than a "loosely bound state". On the other hand, the mass of the 3 Si bound state 
is naturally expected to be higher than the 1 Sq bound state due to the hyperfine-splitting interaction. Indeed, we 
observe that the ground state in the 3 Si channel is much closer to the threshold energy as \ AE\ w M e /25. Although 
the energy level of the ground state is too near the threshold to be simply identified as a bound state, we may expect 
that the 3 iSi ground state is a near-threshold bound state or a loosely bound state. Needless to say, to draw a solid 
conclusion, we need more rigorous signatures of bound-state formation in the 3 Si channel. 
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We employ the diagonalization method to separate the first excited state from the ground state. Fig. [T5] shows 
L-dependence of energies of the ground state and the first excited state in the 3 Si channel for Q = 4. The horizontal 
axis is the spatial size L and the vertical axis is the energy of the ground state (full circles) or the first excited state 
(full diamonds). The horizontal lines represent the threshold energy of the e~e + system together with the 1 standard 
deviation, which is evaluated as twice the measured electron mass. Although it seems that the ground state has no 
appreciable finite-size effect for L larger than 20, the 3 Si ground state lies too close to the threshold energy to be 
assured of bound-state formation. 

As shown in Sec. IIV C[ the distinctive signature of bound states is given by an information of the excited state 
spectra: if a bound state is formed, the lowest (n = 0) scattering state could appear just above the threshold (2M e ), 
but far from the n = 1 energy level of non- interacting two-electron system (2E e {pi)). Indeed, we observe that the 
first excited state appears just above the threshold and its energy rapidly approaches the threshold as spatial size L 
increases. The first excited state can be clearly distinguished from the n = 1 scattering state. Of course, it indicates 
that the ground state should not be the lowest scattering state. Thus, we can conclude: the 3 Si ground state should 
be the S-wave bound state, of which formation clearly induces the sign of the scattering length to change. Therefore, 
the lowest (n = 0) scattering state approaches the threshold from above, the same as the repulsive system in the 
attractive channel. This result shows that our proposal could be quite promising for identifying a near-threshold 
bound state or a loosely bound state such as a hadronic molecular state in a finite box on the lattice. 

B. Bound-state pole condition 

Next, we evaluate the phase a® from the energy level of the ground state through the phase-shift formula J3]) as 
we did in Sec. IIV C 21 All results measured at five different lattice volumes L 3 x 32 are tabulated together with 
results of p 2 and cotoo in Table IXll Indeed, we observe that the phase <jq gradually approaches —45 deg. (— 7r/4) as 
spatial lattice extent L increases. However, <jq is not really close to —45 deg. even at the largest volume (L = 32), in 
comparison to Oo from the smallest volume (L = 16) in the So channel. In this sense, it is hard to judge how large 
of a lattice size is enough to deal with the asymptotic solution of the bound state even in finite volume. Thus, we 
should examine the L-dependence of the specific quantity, cotooj by reference to Eq. [)28p. where the finite volume 
corrections on the bound-state pole condition are theoretically predicted. 

As shown in Fig. [IBJ the values of cotuo are plotted as a function of spatial lattice extent L. Full circles are 
measured value at five different lattice volumes. The solid and dashed curves represent fit results with a single leading 
exponential term and six exponential terms in Eq. (|28[) . The four data points in the region 20 < L < 32 are used for 
those fits. The fitting with the six exponential terms yields a convergent result of 7 as shown in Table IXIII Either 
fit curves in Fig. [To! reproduce all data points except for data at the smallest L. Indeed, the resulting x 2 /d.o.f. is no 
longer reasonable as x 2 /d.o.f. « 3 if the data point at L = 16 is used. Therefore, the ground state at least for L > 20 
can be identified as a bound state without ambiguity. 
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FIG. 15: Energies of the ground state and the first excited state in the 3 Si channel of the e~e + system as a function of spatial 
lattice size for Q — 4. The horizontal line represents the threshold energy determined by 2M e . Full circles and full diamonds 
are measured energies for the ground state and the first excited state respectively. A solid curve with full squares shows twice 
of the single electron energy with non-zero smallest momentum pi — 2n/L. 
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TABLE XI: Summary of the relative momentum squared p 2 , the phase ao and cotoo measured in the 3 Si channel for Q = 4 
at five different lattice volumes L 3 x 32. 





16 


20 


24 


28 


32 


P 2 


-0.0255(26) 


-0.0158(23) 


-0.0167(23) 


-0.0148(26) 


-0.0120(34) 


cro(deg.) 


-56.1(2.9) 


-57.1(4.5) 


-49.2(1.5) 


-47.6(1.2) 


-47.2(1.6) 


cot CTq 


-0.673(73) 


-0.65(11) 


-0.864(45) 


-0.914(38) 


-0.926(53) 


TABLE XII: Fitting results for the bound-state pole condition 


in the 3 Si channel for Q = 4 using Eq. 


(|28|) with variation of 


the number of exponential terms. 










fitting range (L) 


# of exp. terms 




7 


X 2 /d.o.f. 


20-32 


I 




0.1109(60) 


0.24 






2 




0.1218(57) 


0.24 






3 




0.1236(56) 


0.25 






4 




0.1242(55) 


0.25 






5 




0.1251(55) 


0.26 






6 




0.1256(54) 


0.27 






7 




0.1256(54) 


0.27 



VI. SUMMARY AND CONCLUSION 

In this paper, we have discussed signatures of bound-state formation in finite volume via Luscher finite size method. 
Assuming that the phase-shift formula inherits all aspects of the quantum scattering theory, we can propose a novel 
approach to distinguish a "loosely bound state" from the lowest scattering state, which is located below the threshold 
in finite volume in the case of attractive two-particle interaction. According to the quantum scattering theory, the 
S'-wave scattering length is positive (ao > 0) in the attractive channel, if the attraction is not strong enough to give 
rise to a bound state. However, the sign of the scattering length turns out to be opposite (et < 0) once the bound 
state is formed. This fact provides us a distinctive identification of a loosely bound state even in finite volume through 
the observation of the lowest scattering state that is above the threshold. We also reconsider the bound-state pole 
condition in finite volume, based on the phase-shift formula in the Luscher finite size method. We find that the 
bound-state pole condition is fulfilled only in the infinite volume limit, but its modification by finite size corrections 
is exponentially suppressed by the spatial lattice size L. 

To check the above theoretical considerations, we have performed numerical simulations to calculate the positronium 
spectrum in compact scalar QED, where the short-range interaction between an electron and a positron is realized 
in the Higgs phase. We introduce the fictitious Q-charged electron to control the strength of this interparticle force 
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FIG. 16: cot ao in the 3 Si channel for Q = 4 as a function of the spatial lattice size L. The solid (dashed) curve represents a 
fitting result using Eq. (|28|l with only a leading exponential term (six exponential terms). 
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and then can adjust the charge Q to give rise to the S'-wave bound states such as 1 S'o and 3 <Si positronium. Wc 
choose two parameter sets (Q, n) that lead to an unbound e _ e + system (Q = 3) and a bound e~e + system (Q = 4) at 
approximately the same mass of a single electron. We observe the following signatures of the bound-state formation, 
some of which arc related to our theoretical proposals, in our numerical simulations. 

• The lowest scattering state has better overlap with the wall-wall correlator than the point-point correlator. This 
tendency is inverted in the case of the bound state. 

• The sign of the .S'-wave scattering length turns out to be opposite (repulsive-like) even in the attractive channel, 
once the bound state is formed. 

• In the bound system, the phase <7o, which is related to the scattering phase Sq(p) and analytically continued 
into the complex p-plane, is near —45 deg. (— 7r/4) which is associated with the pole condition of the 5*-matrix. 

• The deviation from the pole condition, cotoo = —1, in finite volume is well described by a finite series of 
exponentially convergent terms with respect to the spatial extent L scaled by the binding momentum 7. 

In particular, we regard the second point, the bound-state formation induces the sign of the scattering length to 
be changed, as crucially important for identifying a "loosely bound state" . This is because one can distinguish it 
from the lowest scattering state even in a single simulation at fixed L through determination of whether the second 
lowest energy state appears just above the threshold or near the n = 1 energy level of non-interacting two-particle 
system. We also emphasis that Luscher's phase-shift formula properly reflects one of the most essential features of 
the quantum scattering theory, namely Levinson's theorem. 
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Kinematics of two-particle states on the lattice should be sensitive to choice of the spatial boundary condition. In 
many literatures, this particular point is often discussed and sometimes applied to explore hadronic or non-leptonic 
decay processes [37], IH| or to search exotic hadrons [H, H(| on the lattice. The main point is that the total energy 
of two-particle states, which is roughly estimated by a sum of the energy of non-interacting two particles, depends 
on the spatial size L unless the relative momentum of two particles is zero. Of course, this is because all momenta 
on the lattice are discretized in units of 2-k/L. Here, we have considered the e~e + system. The total energy of 
electron-positron states is approximately estimated by using the naive relativistic dispersion relation in the following. 



zero relative-momentum is not kinematically allowed, so that the lowest energy of two-particle scattering states is 
expected to be very sensitive to the spatial lattice size. In other words, different types of spatial boundary conditions 
(periodic or anti-periodic) exhibit different energy levels of the two-particle scattering states even at the fixed spatial 
size, while a mass of e~e + bound states (positronium states) should be insensitive to the spatial boundary condition 
for the electron fields (48|. For an example, in the n = case, we obtain 



Acknowledgments 



Appendix A: Sensitivity of mass spectra to spatial boundary conditions 




2M e 



(P.B.C. for all spatial directions), 




2E e ( Po ) 



(34) 



2J A/,? + 3 • (I) 2 (A.P.B.C. for all spatial directions), 



which imply an inequality E° e (P.B.C.) < E° e (A.P.B.C). In general, we can expect that ££ e (P.B.C.) < ££ e (A.P.B.C.) 
is always fulfilled. 
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Recently, such sensitivity of spatial boundary condition is often utilized to dis ting uish between two-particle scat- 
tering states and a single-particle state (a bound state or a resonance state) [U [42[. However, there is no rigorous 
test of whether this approach is adequate for such purpose so far. In this subsection, we examine this approach in 
our simulated e~e + system. 

We use the following operators under the anti-periodic spatial boundary condition for electron fields: 

Qp ® = ^^*(x,i) 75 *(x,t), (35) 

X 

*W*) = ^^W(y,f) 75 *(x ) t) e ^-( x -y) ) (36) 

= ^^*(y,<)75*(x,i) e lpi '( x - y ), (37) 

where po = f;(l, 1, 1) and pi = ~(3, 1, 1). The first operator is a simple local-type operator. The second and third 
operators project both the electron and the positron on non-zero lowest momentum (|po| = : ^ L ) and non-zero second 
lowest momentum (|pi| = -^J 2 -), respectively. We can expect that n — and n = 1 scattering states have strong 
overlap with SIm arL d ^Mi j while the bound state has the better overlap with p than £Im an d ^Mi ■ 

Figs. [17] show the effective mass plots of the PP, MqMq and M\M\ correlators in simulations at L — 28 for Q = 3 
(the left panel) and Q = 4 (the right panel) in the 1 S'o channel. There is a similarity between Figs. [5] (P.B.C.) and 
Figs. [17] (A.P.B.C.). Very clear plateaus are given by the MqMq correlator in the Q = 3 case and the PP correlator 
in the Q = 4 case in Figs. |T7l while the same quality shows up for the WW correlator in the Q = 3 case and the 
PP correlator in the Q — 4 case in Figs. [5] In either P.B.C. and A.P.B.C. cases, the PP correlator strongly overlap 
with the Q = 4 ground state, which has already been identified as the bound state in Sec. [IV] Both WW and MqMq 
correlators are expected to have large overlap with the lowest (n — 0) scattering state under each spatial boundary 
condition. A main difference between Figs. (P.B.C.) and Figs. [T7] (A.P.B.C.) is that the PP {MiMi) correlator for 
Q = 3 (Q = 4) in the A.P.B.C. approaches the plateau much faster than the P.B.C. cases. This is simply because 
.©"(A.P.B.C.) is larger than ©"(P.B.C.) and then propagations of non-ground state can die out more quickly in the 
A.P.B.C. case than the P.B.C. case. Indeed, the M\M\ correlator approaches the plateau faster than the MqMq 
correlator in the left panel (Q = 4) of Fig. [T7l since the M\M\ correlator hardly overlaps with the n = scattering 
state as shown in the right panel (Q = 3) of Fig. [171 

We finally employ the diagonalization method to extract the ground states in both Q — 3 and Q = 4 through the 
same procedure described in Sec. IIVB1 Then we compare ground state energies for both Q — 3 and Q = 4 in the 
P.B.C. with those in the A.P.B.C. in Figs.[TS] The left panel (Q = 3), the effective mass of the ground state is clearly 
shifted up in changing from P.B.C. to A.P.B.C, while the plateau of the ground state doesn't change between P.B.C. 
and A.P.B.C. cases in the right panel (Q = 4). An energy shift in the Q — 3 case is consistent with an estimation 
of 2(©° (A.P.B.C.) — ©° (P.B.C.)). We certainly confirm that the scattering state (Q = 3) is sensitive to the spatial 
boundary condition, while the bound state has no dependence of the spatial boundary condition for the electron fields. 




FIG. 17: The effective masses in the 1 So channel on the lattice with L = 28 as a function of the time-slice t in lattice units 
The anti-periodic boundary condition is imposed in the spatial direction. The left (right) panel is for Q — 3 (Q — 4) electron 
fields. Full circles, full squares and full diamonds are obtained from PP, MqMq and M\M\ correlators respectively. 
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FIG. 18: The effective mass plots for the largest eigenvalue X a (t,to = 7), which corresponds to the lowest energy state, on the 
lattice with L = 28 in the So channel. The left (right) panel is for Q — 3 (Q = 4). Full circles (full squares) represent the 
ground states from simulations with periodic (anti-periodic) boundary conditions for electron fields in the spatial direction. 
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Appendix B: Volume dependence of the spectral weight 



The spectral decomposition of the matrix correlator is given by 

a=0 

with the spectral amplitude (u a )j = (0|flj|a)y. A subscript in a ket \a) v stands for finite volume V = L 3 . Remark 
that the finite- volume states \a) v are normalized to unity. The spectral amplitudes are given by solving the following 
equations [H[ 

Y^(w a )i(v* p )i = b af) e E ^l\ (39) 

i 

where (w„)i is an i component of vectors w Q , which are determined through the following generalized eigenvalue 
problem [431 ]: 

G(t)w a = X a {t,t )G(t )w a . (40) 

To solve this eigenvalue equation, we have employed a diagonalization of the transfer matrix M(t,t ) = 
G- 1 ' /2 (t a )G(t)G- 1 / 2 (t ), which provides the same eigenvalues \ a (t,t ) = e - B <*(*-*o) f Eq.®]) 

M(t,t )vL a = \ a (t,t )u a (41) 

with the orthonormal eigenvectors u a = G 1//2 (£o)w Q , if G(t) is an Hermite matrix [43|]. The relative overlap between 
the chosen operator and energy eigenstates (a = 0, 1, 2, ■ • ■) can be determined by the squared normalized amplitudes 
(the normalized spectral weights) 

( ^' = OTF' (42) 

The normalized spectral weights calculated in simulations at L — 28 are tabulated in Table IVl as typical examples. 

Here, we remind that the finite- volume states are normalized to unity, regardless of whether the single particle state 
or the multi-particle state. Suppose the eigenstate a is a single particle state, we simply obtain the correspondence 
between the finite-volume and the infinite-volume states: 

a )oo ( sm gl e particle state), (43) 



where \cx) is normalized as 2E a particles per unit volume. On the other hand, if the eigenstate a is a two-particle 
state, the correspondence factor between la)^ and \a) v should depend on dynamics between two particles. Such 
corrected factor is explicitly derived by Lellouch and Liischer (denoted in the following by LL) to determine the 
physical K — > tttt amplitude from the finite volume calculation [351 ]. Here, we consider the LL-factor only in the 
non-interacting case, where the scattering phase shift between two particles is taken to be zero, for a simplicity, and 
then obtain the following correspondence between la)^ and \a) v for S'-wave two-particle states [35l [36j| : 

a )oo (two-particle state), (44) 



E a V 



which indicates that observed spectral amplitudes (v a )i for the local-type operator (i — P) are proportional to 1/VV 
for the single particle state and proportional to l/V for the S'-wave two-particle state, since the physical spectral 
amplitude (0|f2p|a) in the case of the local- type operator should not depend on the size of the spatial volume 
V = L 3 . 

Let us consider the volume dependence of the spectral amplitude of the ground state with the local-type operator 
Hp. In Fig. mi we plot the finite-volume spectral weight |(0|Slp|a = 0) v \ 2 scaled by V 2 for Q = 3 and by V for Q = 4 
as a function of spatial lattice size L. Recall that Q = 3 is the unbound system, while Q = 4 is the bound system. 
No appreciable L-dependence is observed in either cases. This indicates that the finite-volume spectral weight for the 
local-type operator has a specific volume dependence according to whether the single particle state or the two-particle 
state. In other words, each contribution from two-particle states (scattering states) relative to the single particle state 
is suppressed by a inverse of the volume factor, 1/L 3 in the PP correlator. 
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FIG. 19: The volume dependence of the spectral weight of the So ground state with the local-type operator fip. We plot 
\(0\Q P \a = 0)| 2 scaled by V 2 = L 6 for Q = 3 (full squares) and by V = L 3 for Q = 4 (full circles) as a function of spatial 
lattice size L. Remind that Q = 3 (Q — 4) is the unbound system (the bound system), where the n — scattering state (the 
bound state) is the ground state. There is no appreciable L-dependence in either cases. 
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